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ABSTRACT 


Several techniques for modeling the disturbances to a space- 
craft's attitude caused by moving crew anoibers are presented. These 
disturbances can be the largest mommts acting on a manned space- 
craft, and knowledge of their effect is important in the sizing, 
design, and analysis/si.;>ulation of spacecraft attitude control 
systems. 

The modeling techniques are identified as two principal types: 
deteministic and stochastic. Ibree techniques of each type are 
pzosente-.. The deterministic models include point^sass motion 
derivatives and a discussion on dynamic models of moving crew members. 
Ihe stochastic t -chnlques are hi^ili£^ted by a Fourier transform 
method and the representation of long-term crew disturbance activities 
as o\itputs from appropriately designed filters. A z-transform technique 
is developed to obtain a difference-equation form of stochastic models 
for use on digital computers. 

An appendix derives spacecraft equations of notion idiich can be 
used with many of the models discussed. 
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INTRODUCTION 


Advances in manned space flight have led to the participation 
of man as a usef\il part of missions involving the carrying out of 
sophisticated experiment programs. He can now contribute to the 
success of missions requiring maximum use of his capabilities as an 
observer, operator, repairman, or subject. His presence in space, 
however, is not without its penalties (which should be more than 
offset by its advantages). Life support systems must be provided, 
at significant costs in weight and spacecraft complexity. Other 
spacecraft systems must be made highly reliable, and redundancy 
provided, in order to insure that man will return from his missions. 
Again, cost in weight and complexity are incurred. Finally, man by 
his presence imposes a disturbence environment on the spacecraft 
which could hinder the carrying out of certain experiment tasks. 

The expected dynamic environi^ent is of primary importance in the de~ 
sign and analysis of any spacecraft stabilizaion and control system. 

For an orbital spacecraft, external distiirbances con be caused by 
aerodynamic forces, gravity gradient torques, or meteoroid impacts, 
for example. The internally-generated disturbances are due to, ?or 
example, rotating component imbeilance, mass (cargo or equipment) shifts 
or transfers, or as mentioned above, by motions eind actions of the crew 
members* I 

Crew motions can produce some of the largest disturbances act- 
ing on a manned spacecraft. These motions can include restrained 
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activities (where an astroimut is performing tasks at a console, for 
example), or translation, where the astronaut may be movirj/ 
point to another within the spacecraft. As motion occurs ^uile 
a crew member is in contact with the spacecraft, ‘disturbance mcments 
are generated. This dynamic interaction of man with his spacecraft 
requires that the attitude control system be designed to take into 
account the disturbance torques arising from the man’s motion within 
the vehicle. 

Early work on "modeling crew motion distxirbances was, in general, 
confined to simple representations of motions within a spacecraft 
(see refs 1-5) as instantaneous or (usually) constant velocity 
translations. This type of modeling is best suited to rigid-body 
spacecraft considerations, as there is no discrete frequency content 
expressed in the movements. Also, no treatment of the restrained 
crew activities is afforded. 

Subsequent investigations have resiilted in the modeling of crew 
disturbance forces and moments as stochastic processes (see refs 6-8) 
with finite band-width, so that frequency domain control system res- 
ponse problems can be explored. 

More recently, determininstic modeling of astronaut body motions, 
and the resultant disturbance forces and moments, has been carried 
out by Kurzhals and Reynolds (ref. 9* Appendix B) using models of the 
human body developed by Hanavan (ref. 10) and modified by Woolley 
(ref. 9» Appendix A)* 
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This thesis focuses on some of the modeling techniques devel- 
oped to date, and indicates some logical extensions to some of 

these. In particular, since much analysis and simulation work is 
carried out by high-speed digital computers, the modeling of sto- 
chastic processes in this domain is deemed important j development 
of difference-equation descriptions of a stochastic crew motion 
model is presented as a potentially powerful technique for computer 
applications. 

TYPES OF CREW MOTION DISTURBANCE MODELS 
Two types of astronaut crew disturbance models will be considered 
herein: deterministic and stocheistic. The deterministic models are 

completely specified as functions of time, eind may rsuige from 
representations of simple point-mass motions to sophisticated dyncmic 
human body models. The stochastic models will generally be rep. ed 

by their statistics (e.g. statistical moments) and frequency content 
or bandwidth. 

Deterministic models will be considered first. 

DETERMINISTIC MODELS 

Crew Motions Represented by Moving Point Masses 

As pointed out in the Introduction, one method of representing 
crew motion disturbances is to treat one or more crew members as 
point masses, and allow this mass (masses) to translate within the 
spacecraft. The simplest means of implementing this method in a 
dynamic analysis is to consider only a single mass (which could 



represent any munber of crew) and exaiuine the effects of its 
motion about the spacecraft. Appendix A describes the derivation of 
this approach. To include the effects of separate masses in motion 
(either simultaneously or in series), modification of the basic 
equations given in the Appt.»dix must be made. For i masses, i=l,... 
total number of moving masses to be considered (n), define 


_ m. (m - m. ) , . 


m 


Equations (A36) of the Appendix would then be written as 
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Moments and products of inertia in the above equations are given by 


I » I + y' Q. /y? ♦ z?\ 
X x,o xj 
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ix * = 2 f A » *i‘i) 
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For these equations, the position tine histories for the n sepenrate 


masses must he specified. 
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Crev Motions Represented t>y a Dynamic Model of Man 


In the previoiis section » where a crewman was treated simply as a 
point mass, no condideration was given to the fact that man possesses 
local moments of inertia, cuid articxilated body segments which, by 
their motion, result in the generation of local forces and moments. 

A more sophisticated means of modeling crew motions is through 
the use of mathematical equations describing the motion of these 
airticulated body members. A simple example will be used to 
demonstrate this concept* 


If it is assumed that an astronaut is in contact with the spacecraft 
at his feet (see figure l), then any forces and moments he generates 
by moving various parts of his body will be transmitted to the 
spacecraft through this point of contact. Assume that the man depic- 
ted in figure 1 is rigid, with the capability of moving his right 
am in the yz-plane. Considering the right arm to be rigid, and 
allowing it to move through the angle 6, force and moment components 
required to produce the movement are given by 



m y 

= a cm 


( 5 ) 



m z 
a cm 


( 6 ) 



-I e 

a 


s 


(7) 
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Figure 1 - Schematic of dynamic model of man in 
contact with spacecraft. 


where 


Y = £ sin 8 (8) 

cm 

Zc^ *-l(l - cos 6) (9) 

• • •• 

The quantities and Z can be expressed as functions of 6 

cm cm 

and its derivatives as 

.» •» • 

^ “ 6*sin6) (10) 

*» M • 

Zgjjj = - i(6 sin 0+ 6*cos 6) (11) 

It can be seen that expressing the angle 6 as a function of 
time allows the computation of the forces and moments required to 
produce this pendulum-like motion, and hence (since the astronaut is 
in contact with the spacecraft) the reaction forces and moments 
acting on the spacecraft. 


F = rn s, {0*sin6 - 0 cos 6) (l2) 

J S ^ 

» •* 

= m i (0^cos6 + 6 sin6) (13) 


M.. 




(iM 


Similar expressions can be obtained for movements of other body 
members. For move complex motions ("double-pendulum" arm movements, 
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for exat-ple), transformations to a single reference axis system are 
req,uired. Appendices A and B of reference 9 present the development 
of equations for describing the notion of a nine-segment model of 
man. Inputs required for the use of this model are time histories 
of Euler angles describing the orientations of the nine body segsents; 
from these inputs applied forces and moments are obtained. 

The dynamic model developed in reference 9 may also be used for 
an astronaut moving vith respect to a spacecraift (as in the perform- 
ance of translational "soaring" maneuvers, for example). For this 
type of application, equations describing the motion of his center- 
of -mass vith respect to the spacecraft are required, in addition to 
the Euler angle histories for the body segments. 

One method of defining the motions of the model segments is to 
record the segment motions of an actual subject performing typical 
activities. This can be done by talcing motion pictures and 
extracting the required information from the resulting film 
data (see, for example, reference 11.) More directly, use of 
an exoskeletal device such as described in references 12 and 13, 
can be made. Potentiometers at the Joints between principal body 
segments permit the recording of relative segment motions on magnetic 
tape or punched cards. These motion histories, with appropriate 
transformations to reference co-ordination systems, can then be used 
to generate applied force and moment histories for input to space- 
craft control system analyses. 
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Crev Motions Represented by Their 
Force and Moment Histories 


For control systems analyses and simulations » it is generally 
most convenient to represent crew disturbance inputs as forces and 
moments along or about the spacecrafts axes. This can be accoiqplished 
(for restrained activities) by attaching a subject to force and moment 
measuring transducers, and recording the outputs of the transducers 
as the subject performs typical tasks. Load cells or a strain gage 
balance are devices which can be used to obtain the force and moment 
histories. A load cell arrangement which has been used successfully 
in both laboratory and fli^t experiments is shown in Figure 2. This 
array features essentially equal sensitivity for both forces and 
moments with respect to all axes. The resolved forces and moments 
are given by the following equations (taking load cell compression 
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Figure 2 - Load cell array for experimental force 
and moment measurement. 
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In simulating (on earth) the forces and moments caused by a 
moving crew member, the introduction of certain moment contributions 
from a subject's center of gravity shift may occur. Consider figure 
3 representing a subject whose disturbance forces and moments (from 
typical limb movements) are being measured. 

The applied forces and moments are measured with respect to the 
point 0, which is on the z* axis. When the subject moves a particular 
body member so as to cause his total c.g. to shift, a moment about the 
point 0 will be produced by any motion of the c.g. other than along the 
z axis. This moment will be equal to the static weight of the subject 
times the distance through which .ne c.g. travels in the x-y plane. For 
example, if the subject raises his arm over his head, his total c.g. will 
move along the (-x) axis, creating a gravity moment about the y-axis equal 
to the -X ' W. Similarly, if he holds his arm at a 135® angle with re- 
spect to his torso (in the x-y plane of figure 3), moments about the y- 
and x-axes will be produced. 

There are several methods by which these gravity moments can 
be mechanically eliminated from the data. These methods consist of 
supporting various body members to exactly balance out the gravity 
force. Such arrangements as suspension slings and underwater neutral 
buoyancy weighting are two examples. However, it is not always prac- 
tical to provide these mechanical suppoirts ( for reasons of available 
space or cost, for example). It would thus be desirable to analy- 
tically remove these gravity-produced moments from resultant 
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Figure 3 - Schematic of subject force and 
moment measurement. 
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disturbance data, keeping in mind, of course, that the presence of 
gravity cannot be eliminated and will therefore inherently affect 
the way in which simulated motions are performed. One method of 
mathematically subtracting out gravity moments consists of 
instrumenting a subject with an exoskeletal device which measures 
all limb rotations (as described in a previous section). These 
rotations can be used, along with mass and center of gravity data 
for each limb, to compute the amount of overall c.g. travel due 
to the articulation of one or more limbs. This method has been 
used (ref. 13) with results ranging from fair to good. In no 
cases were perfect fits obtained, however, since the limb motion 
measuring device covild not be fitted to measure every conceivable 
body member motion. 

Another method of reducing applied moment and force data to 
eliminate gravity contributions is described in the following 
analysis. For this derivation it is necessaiy to assume an axis 
system orientation; the orientation shown in figure 3 will be select- 
ed, that is, gravity forces are exerted in the +z direction. 

As the subject performs various motions, forces and moments are 
produced at the point 0 in the force measuring unit (which covild 
consist of load cells or a strain-gage balance). These forces and 
moments are measured, and are as follows: 
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Since the subject is in contact with the force measureing unit, 
there is exerted on the subject equal and opposite reaction forces 
which produce accelerations of the body c.g. 




The mass of the subject is simply equal to his weight divided 
by the acceleration of gravity constant: 


m 

s 


W 

s 

g 


(18) 
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To cos$)ute the c.g. travel then, a double integration of the forces 
may be used. 


XI 


P dt dt + jc t + X 
X o o ^ 

s c.g. c.g. 


6 * 


c.g. 


•g* 


Since the subject will normally b«;gin all simulation sequences 
from some calibiated rest position, the initial velocities and 
desplacements of the c.g. can be taken as 0. Thus 

( t) = _g_/yP„ dt dt = - __g_ /•/•F„ dt dt 




m 




C.g. 


( t ) = g F dt dt = - 




tff 


F dt dt 
’F dt dt 


The correction for c.g. shift would then be written as follows 


M = M 

X X 
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measured m V JJ 

M *M-x W*M+gW /•/•?( 

^actual W JJ *ni 


dt it 


dt dt 


M * M 

z z 

actual m 


>(19) 


(20) 


( 21 ) 
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Remember that these corrections are for the parti* ar axis system 
orientations assumed. It is seen that the mat if .ical correction 
is independent of subject weight (which appears implicitly in the 
variables) and dependent on the j'lective gravity acceleration 
constant at the location of the simulation. The correction equations 
could be used at any level of apparent gravity by proper adjustment 
of g. Thus, in a zero-g orbital, environment, where gravity 
acceleration is exactly counterbalanced by centrifugal force, the 
effective g would be zero, and measured moments would be those 
actu€dly produced by the crev motion. 
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STOCHASTIC CREW MOTION MODELS 


In many areas of spacecraft control system analysis, there 
is either no requirement for deterministic representations of crew 
disturbance forces and moments, or their use might be unwieldy. 

Instead, the disturbance forces and moments can be represented by 
other characteristics, such as frequency content or statistical 
mcHnents. The so-called stochastic crew motion disturbance models 
can be used in the time domain for control and pointing system 
simulation analyses, or in the frequency dcmiain for systOL sizing 
or structural analysis problems. 

Fourier Series Representation of Force and Moment Histories 
If time histories of simulated crew motion distxirbance forces 
and moments are available (for selected activities), then a Fourier 
series representation of each force and moment history can sometimes 
conveniently be obtained. 

For a force (or moment) history, F(t) over a time interval t = 0 to 

t » t^ there can be derived Fourier coefficients A and B 
f n n 

(n * 0, 1, 2, ...) to obtain the Fourier series 
A * 

f(ti = — V* [A cos n u t ♦ B sin nw t] (for 0 < t < t.) 

^ 2 + / ^ n on o i 

n=l 


( 22 ) 



where 


■ t /:- 


(t) cos dt 


and 


B 


■e/' 


(t) sin Do> t dt 
o 


for n = 0, 1* 2, 


for n -= 0, 1* 2^ 


.« (23) 


.« (2U) 


•Ihe quantity (*)^ is the fundamental frequency of the series and 
n is the harmonic of this fundamental frequency. In actual practice, 
n would have seme finite limit, ^ich vo\xld depend, in part, on the 
error limits desired between the Fourier series and the original data. 
The fundamental frequency is taken to be 

- 2 TT 

•"o ■ t^- (25) 


A set of six Fourier series can thus be obtained for the force 
and moment time histories representative of a given period of crew 
activity. A constraint which should be noted is 


/* V 

(t) dt = ’ "I 


dt = 0 


( 26 ) 


o o o 

A non-zero integral wo^old tend to produce a net change in momentum 


of the man-spacecraft system, if the force and moment models were 
applied to an attitude control system simulation for example. 
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In addition, it should be assumed that the activity would be 
initiated from (and return to) some condition of rest, i.e. net forces 
and moments equal to zero. This condition requires the additional 
constraint that 


A 

O 

2 


♦ 



A 


n 


0 


n=l for tdiatever H chosen. 

In practice, the number of Fourier coefficients required (i.e. the 
veiLue of N) may be quite large to give a good ^proxiaation to 
measured data, unless the measured wave form is periodic with a 
relatively saaJ.1 bandwidth. For those activities whose force and 
BKMBent histories are such that the foregoing Fourier series models 
are not readily derivable, more formal Fourier transform techniques 
may be employed. 


Use of Fourier Transforms in Crew Motion Modeling 
A time history (of, say, a force or moment component of a crew 
motion disturbance), P(t) has a Fourier transform 


f((u) e 



( 27 ) 


and since F(t) can be assumed to be zero for all time except U ^ t ^ t^. 


f(u>) = 



(28) 
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This equation is essentially equivalent to the tvo eqm^tions 


used to compute the A^*s and in the proceeding section » stn^e 


-iu)t 

e 


cos (ut • i sin (Dt 


< 29 ) 


It is used vhen the force or aoaeat conponent possesses a contxauous 
frequency spectrum* rather than a few discrete frequencies as is the 
case in the preceeding Fourier series disciission. In practical 
situations involving crew motion disturbance forces and mnents* 
frequency content is limited to a range, say 

|(i)| < |u)j,| , idiere <o^ may be arbitrary. 


Thus f(b>) would be essentially zero for all Nl ^ and the time 

functions, F(t), can be expressed eis (from ref. 6). 


F(t) 


1 

2it 


/ 




f(<i))e^*^du> 


( 30 ) 




These relations, f(u) and F(t), constitute a Fourier transfoim pair, 
and apply to continuous wave f'^rms, and thus ^lave somewhat limited 
use (e.g. , in analog cimputer simulations or analyses). 

For most applications Involving control system analyses (at 
least in recent years), the high-speed digital computer is utilized. 
This in turn requires tnat discrete versions of the Fourier transform 
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be used, and, of course, all wave foras or functions oust be sampled 
and not continuous. 

The discrete Fourier transform pair associated with f((t>) and F(t) 
respectively are 


and 


H 


f{j) • S P(k)e'' 




(j » 1,2,. ..F) 


(31) 


F(k) « ^ f(j)e 


J=1 


(k * 1, 2, ...N) 


(32) 


where H is the number of samples in the wavefora. F^k) and f(j) 
are generally complex sequences defining functions in the time and 
frequency domains respectively. The index k refers to a sanqple 
period nusdier from the time history, and the index J represents 
a hanaonic of the basic frequency in the frequency function. The 
time and frequency corresponding to these indices can be found from 
the relations 

tj^ = (k-1) At I where At = | 


fj * (j-1) Af = (j-1) ^ |wbere Au * 

A note of caution is required. The sequences f(j) and P(k) are 
periodic over H because of the exponential in equations (31) and 
(32). la obtaining the response of a control system, for exanqile. 
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this periodicity becomes of concern when the force or moment input 
is convolved with the transfer function of the system under study. 
Consider the following application of crew disturbance modeling. 

A tine history representing disturbance forces and moments will 
be a real function; hence, the sequence F(k) will be real, and 
the sequence f(j) of the discrete Fourier transform will be complex 
and exhibit the properties that 

(1) Re f(j) = Re f(N+2-j) (j*2, 3, ...|) 

and 

(2) Im f(j) = -Im f(N+2-j) (j=2, 3. ...|) 

That is 9 the real part of f(j) is symmetric about the frequency 
0)^ = — idiich is referred to as the "folding frequency," and the 
imaginary part of f(J) is antisymmetric about u^. Since the 
sequence f(j) is periodic. Re f(j) can be considered an even 
function, and Im f(j) an odd function. The Foxirier coefficients 

II 

f(^ 2)...f(N) can be thought of as the negative frequency harmonics 

for frequencies between - (u)_ - u ) and -ui . 

r o o 

For the time history being used, enough sasq>les, N, must be available 

to insure that all expected frequency content of the input is at 

frequencies lower than the folding frequency . Also, the 

^ ^f 

transfer function of the control system being studied, H(s), (s * iu)), 
mvist be adequately described between u» = and u = If the input, 

F(t) is not periodic over the interval t = o to t = t^, then eui 
Incorrect system response r(t) will be obtained from convolution 
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of the input with the transfer function. To obtain the correct 
response an additional N terms of zero magnitude should be added to 
the input sequence to obtain a sequence of 2N terms. The discrete 
Fourier transform* the inverse transform* and hence the convolution 
involving this 2H > term sequence will then be periodic over 2H 
points* and an N~point convolution of ^(h} and H(s) will yield 
the aperiodic convolution of two N-point sequences. The transfer 
function* H(m)* is defined as 

H(m) = H [i (m-1)] m * 1, 2*...H+1 

over the first (N+1) points, and 

H (N + 1 + j) = conj H(N + 1 - J) J » 1, 2*...H-1 

over the remaining (N-1) points* for a total of 2M points in the 
sequmce. Convolving li(m) and f(m) over m* 2N i>oints gives the 
frequency response array 

R(m) H(m) f(m) m » 1, 2*...2N 

The inverse discrete Fourier transform given by equation 32 can be 
applied to R(m) and the sampled output response* f(t)* can be 
obtained from 

r (m) ■ Re R(m) m » 1* 2*...N 

by noting that 

t • (m-1) dt = (m - 1) ^ 

N 

The system response to the crew motion disturbance input over the 
interval 0 ^ t ^ t^ has then been obtained. 
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Crev Motion Disturbance ApproxlEation by Filtered Randoa Hoise 


A great majority of astronaut crev motions have been, and 
undoubtedly will continue to be, associated wi(.h long-term, 
essentially continuous, activities such as experiment console monito- 
ring, EVA preparations, meal preparations, or exercising. These 
activities, or ccsnbinations of activities, are composed of more or 
less random body movements, resultant forces and moments thus 

appear somewhat random in nature, and resemble stochastic processes. 
This section will discuss the modeling of crew motion disturbance 
forces and moments (for continuous activities) by treating them 
as stochastic processes. 

A stochastic process can be defined as a family of functions of 

time, y^(t), the realizacion of which depends on eua "outcome” (which 

can occur with some probability p) of some "experiment" (ref. lU). 

A relevant example mi^t be the force applied edong the y-axis by an 

astronaut perfoming a task such as preparation of a meal or taking 

a shower. The first time the task is performed (taking some length 

of time, T) the applied force history could be represented as F (t). 

^1 

The second performance would probably involve sli^tly different 

movements, and the force history F (t) would be different, in a 

^2 

deterministic sense, from F (t). Similarly, F (t), F (t), etc. 

Fi 

would be expected to be different from each other. For a single 
performeuxce then (a given i) P (t) denotes a single time function. 

y 

For a specific time t. for any performance F is a random variable, 

u y 
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dependent upon the performance i. Finally, for a given performance 

at a given time, F (t ) is a single number. 

j 

Since a stochastic process can be thought of as a time series 
of random variables then the statistics of the random variable (such 
as mean, variance, autocorrelation, etc) can be used to characterize 
the process. 

If the random process under consideration is stationary (i.e. , 

its statistics, such as mean, variance,etc. ,are not affected by 

a shift in the time axis), the statistics will not be affected by a 

shift in the origin of the time axis. 

Further, if the statistics of a reuidom process can be determined 

from a single function (in the present exanq>le F (t) ) of the 

^i 

process, then the process is said to be ergodic, and the various 
statistical parameters are expressible as time averages. For a 
stationary random process (again use F (t) ) the power spectrum 

y 

is the Fourier transform of its autocorrelation R(t) given as 

S(u)) » R (t) dT (33) 

00 

>diere (from ref. 15) 


R (T) 


/ T/2 

Fy(t) . Fy (t ♦ t) dt 
T SB - T/2 


(3V) 


and T is a time "leg" between two sanqples of the time function 
F (t). It should be noted that the eq'iation defining R(t) assumes 

y 
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that the mean of the function 


P (t) is zero, or 


1 

T 




o 


dt « 0 


a constraint idiich is essentially the same as that discussed in a 
previous section (see eqn 26). 

Once the pover spectrum has been obtained for a given force or 
moment history, appropriate mathoDatical filters can be developed 
so that a filter output, vhen excited by scxne standard input (random 
or "white" noise is convenient to use), possesses the same pover 
spectrum characteristics as the force or moment history being modeled. 
In equation form 

" ^filter * ^Input (35) 

using the example of the force F^. 

The next section presents a developme..c, using z-transform 
techniques, of an extension of this model. 
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A z-TRANSFORM TECHNIQUE FOR CREW MOTION MCOELING 

As previously noted, many crew motion disturbance modeling techniques 
will find use in digital computer simulations and analyses of spacecraft 
stabilization and control systems. Digital computer representations of 
real-world differential equations used to describe physical systems makes 
use of difference equations, in contrast to the use of integrator ele- 
ments in an analog computer. Formulatxon of these difference equations 
for use on the digital computer can, in many cases, be most easily accom- 
plished by eiq>loying z-transform methods. The purpose of this section is 
to develop the methodology for generation of difference -equation repre- 
sentations of crew motion disturbance models, through use of z-transform 
techniques . 

If f(t) is a continuous function (in the time domain), then the 
z-transform of the sequence f(n At), obtained by sampling f(t) every 
At seconds is defined as 


Z[f(n At)] = F(z) ■ £ f(n At)z”*^ (36) 

n*0 

where z is a complex variable and 



(37) 


where p is the radius of convergence of the series f(n At)z ” (ref. l6)). 
It is noted that the z-transform can be related to the Laplace transform 
by 
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8 At 


(38) 


z ® e 

where s is the Laplace variable. 

Reconsider the filtered random noise representation of crew 
motion disturbances which was outlined in an earlier section of this 
paper. It is recalled that Fourier transforms were utilized to repre- 
sent the filter structure, with power spectra computed in the frequency 
domain. It has been noted (refs. 17 » 18, 19) that Laplace (or Fourier) 
transforms of sampled functions result in infinite series, %rhich are 
nonalgebraic (containing factors of the form e® ^^), making handling 
of sampled-data, or discrete, systems more involved and often unwieldy. 
Use of the z-transform in describing the transfer function of the model 
filters enables the derivation of difference equations which are readily 
programable on the digital ccmiputer. The following discussion develops 
the applicable theoty for obtaining the requisite difference equations. 

It should be noted that this development is most amenable for use with 
activities of an essentially random nature, whose force and moment time 
histories would be expected to possess relatively smooth power spectra. 
Short-term, discrete, activities would be more suited to modeling by 
other, dete^'ilnistlc techniques. 

Assume a set of N data points as a result of obtaining a force or 
moment time history (M) from some laboratory or experimental crew motion 
disturbance measurement activity. . The sampling period, which will be 
assumed to be constant, is at. A third assumption which will be made is 
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that the data, represented hy the N data points, has zero mean. Fol- 
lowing Blackman and Tukey (refs. 13 and 2u), the correlation function is 
first deterained as 


, T/2 

R(t) = lim ^ / M(t)M(t + T)dt (39) 

-T/2 

T -► o> ' 

which, for the discrete case, can be written 

1 ” 

R(t) » lim ^ E h(k)M(k + t) (UO) 

N-*-« -^k = 0 

where T is commonly referred to as the "lag" and is an integer multi- 
ple of At (t = n At). 

The one-sided power spectral density (PSD) is <;iven as the discrete 
Fotirier transform of the autocorrelation 

00 

S(u) = 2 At Z R( nAx)e"^‘^^‘^ 0 £ w < « (Ul) 

n * -*“ 

The PSD, S(u), is a real, pcsitive function of u>. This is the 
function to be modeled. 

The next step is to assume a transfer fiuictlon, H(s), which will 
be excited by a random, or "white," noise input to yield an output with 
the power spectrxm characteristics of the time history, m. H( s) is 
an "open-loop" transfer function, and the output PSD as a function of u 
can be determined in a conventional manner as 
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Output ■ 


input 


(U2) 


where H(iw) ■ H(s)| _ , 

8 * iW 

It is now necessary to define a fom for H(s) which will give '^he 
appropriate output PSD. In the frequency domain, this filter mu''*’ yield 
a zero output at zero frequency, in order that a resultant force or 
m^ent history not exhibit a bias (which would result in addition of net 
momentum to the spacecraft dynamic system). Accordingly, an arbitrary 
transfer function is assumed; following Hendricks and Johnson (ref. 7) 
this function might be specified as 

H(s) * 5 ski (43) 

[[sf - 2b^(s) + c^^l [{sf - 2 b 2 (s) ♦ c/] 

If it is asBianed that a random (Gau .^e input with an rms rilue 

of unity, the model PSD as a function of w is (f.ron equations (U2) and (1*3) ) 

SM 5 ^ 5 5- • 1 (W) 

output 1 + U) ♦ Bg w + w + 0) 

where A and B^(i * 1 - U) are functions of the filter parameters 
a, bi, bp. Cl, Cg. A method of curve- fitting must be used to fit the 
model PSD (equation (U4)) to the measured PS? (equation (4l)). A form 
of least- squares curve-fit algorithm could be utilized to determine the 
factors A and B^, and hence the filter parameters. 
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Once the transfer function has been defined in the frequency domain, 
derivation of an equivalent digital filter is pursured. It should be 
noted that the filter being developed (for use on digital computers) will 
essentially receive input from a digital sampler, whose output can be 
considered a series of impulses. In order to more accxirately represent 
the continuous input to tie continuous filter being modeled, tfc' iero- 
order hold is ccmnonly used (ref. 21). The mathesiatical foranilation for 
this "element" is given by 


».ch<=> 



(l>5) 


The factor — in H , (s) must be included with the filter transfer func- 
s zoh 

tion H(s), and use is made of equation (36) to obtain the digital filter: 


H(z) = Z[H.(s)] 

Z 1 


(U6) 


where 

1 s 

For the present problem, H^(s) becomes 


H (s) = — ^ s- (U7) 

(s - 2bj^s + )(s - 2 b 2 S + ) 


The transfer function H^(s) given in equatio.i (U7) posses poles 
(which will be complex for most cases) at. 
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To obtain the digital filter, H(z), corresponding to the model filter 


(vith zero-order hold) H^(s), Cauchy’s residue theorem (and eqn. 38) can 


be used (ref. 19 ) 


2 ^ 1 1777 


(p) 

At(s - p] 



M 

P 

I 


(p - 




where H = number of poles of H. (s) (in this case N “ U) and p. is the 
p i P * 

value of the kth pole. The path of integration in the z-plane must in- 
clude all poles p . From equation (W), a ratio of polynomials in (z 

K 

can be obtained:^ 


^It would have been as easy to obtain a ratio of polynomials in z; 
however, as will be seen, the given ratio has certain advantages in 
subsequent developments. 

3U 



H(x) 



-1 -2 -1 
Oq X ♦ oig X 0|^x 

6q ♦ x”^ ♦ Bg +...•♦• B^*“* 


(1»9) 


where the a*s end B*s are functions of a, b^, bgt c^* 

The transfer function i^Ewt-output relation can now be applied to 
determine the difference equations to be used on the digital conputer. 
It is first noted from x-transfom theory (refs. l6 and l8) that an 
inverse transform may be expressed as 


»-l 


X* G(x) 


E g(n &t m &t)6(t - n Lt) 
n^O 


(50) 


for the sanple sequence where t >, 0, g(n At -t* m dt) is a discrete 
ssBQile from some time history g(t) whose z-transform is G(z), and 
£(t - n dt) is an impulse function. Equation (U9) is rewritten as 


(6q Bj^ x~^ +...+ Bjij z”“)F(z) = (Oq + z“^ +...♦ aj^x“*’)G(z) 

(51) 

The inverse z-transform of equation (51) is taken term by term to yield 

B^ £ f(n dt)fi(t - n dt) + B, E f f(n - l)dtl 6(t - n dt) 
n«0 n*0 *• ■* 


♦• • •+ B 

E f f(n - m)d^ 6(t - n dt) 

Ja 

n=0 •- ^ 

€0 

* r T 

« Gq E 

g(n dt)6(t - n dt) + a. E gj(n - l)dtj6{t - n dt) 

° n=0 

1 n«0 *- J 


00 r n 


E gj(n - l)dtj 6(t - n dt) (52) 

r,=n •- J 
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Ihe coefficients (B^ f[(n - i)At] and a. g[(n - i)At] of the 

inpulses 6(t - n At) at the sane sampling instant, n At, can he 
equated to get 

B« f(n At) ♦ B, f((n - l)Atl B f[(n - m)At] 

Ox 01 

* g(r. At) + g[(n - l)At] ♦...■*• g[(n - l)At] 

( 53 ) 


which is the difference-equation description of the process obtained by 
passing the "white" noise signal throxigh the filter H. Rewriting 
equation (53) as 


f(n At) = g(n At) + ^ g[(n - l)At] +. 





g[(n - A)At] - f[{n - l)At] 


-^f[(n - 2)At)] 


e 

- ^ f[(n - ai)At] 
^0 


(5h) 


it Is seen that the output f(t), t = nAt, is a function of the backward 
differences of input and output or, in other words, the output at a given 
sasqile instant is dependent on both output and input at earlier sample 
instants. The quantities g(nAt) will be values from the "white" noise 
sequence, and as such, will be merely a set of normally distributed random 
numbers with a mean of zero and a veuriance of unity. 

Implonentation of equation (5^) on a digital computer is straightfor- 
ward. Since the difference equation was developed for t > 0, the values of 
input and output prior to t = 0 (i.e., g(-At), g(-HAt) ,. . . ,f (-At) , f(-tAt), 
etc.) can be taken as zero. 
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Example 


As an illustrative example of the utility of the difference-equation 
description of a filter output, consider a relatively simple filter 
structure to be given by 


' (s+2)(s+l<) 


(55) 


Following equations (U5) and (U6), to include the requisite zero-order 
hold, and equation (U8), H(z) can be obtained as 


H(z) 


z-1 ( 1 r z 

z I 2 I -2At 

y ^z-e 



(56) 


Assuming At - .005 sec., H(z) becomes (after simplification) 

H(2) = » 00^92558025 z~^ - .00U92$$8025 

I-I.9702U8507 + . 9701*1*55927 z‘^ 


(57) 


The difference-equation relation between input, x., and output, x , with 

10 

being given as a sampled input, can be written 

X (nAt) = .001*92558025 x.((n-l) At]-. 001*92558025 x.((n-2)At] 

° ^ . (58) 

+ I.9702U8507 X [(n-l)At] - .970UU55927 X [(n-2)At] 
o o 

To ascertain the accuracy of equation (58) in describing a system, 
output response to a unit step function^ through the filter was 
determined using a small desk-top computer system. Figure 1* shows 
a plot of the equivalent analog response (solid curve) with the 

^The unit step function is defined (ref. 17) as having the value 1 for 
positive time, t, 0 for negative t, and 1/2 for t = 0. 
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Figure Response of exajaple filter to unit step input 




59 


Figure 5.- Plot of difference between analog and difference-equation results. 



difference equation output for comparison. As can be noted from 
figure the difference-equation result is essentially eqiial to the 
analog curve after O.U seconds. Figure 5 gives the difference (analog 
minus difference-equation) between the two results for the same 2.5 
second period shown in figure U. Agreement of the difference-equation 
description with the euialog result is within two percent from O.U 
seconds on. 

Next, response of the filter (equation (5T) ) to a randan noise input 
will be examined for its implementation characteristics. 

The technique used to digitally generate pseudo-random numbers for 
input to this problem consisted of first obtaining pseudo-random numbers 
uniformly distributed on the interval (O, l) through the relation 


Fractional part of [(it + u^^ 


(59) 


Normally-distributed pseudo-random numbers , N^ , with mean 0 and unity 
variance are then generated from (see ref. 22 ) 


Nj^ = (-2 In cos (2 it 

■ (-2 In sin (2ir 


(60) 


Uo 




X 


o 


111 


Response 


where and imist he independent uniform variables. This 

condition necessitates the use of two sequences of the form given by 
equation (59). with separate initializing uniform random numbers. 

The generated are used as step inputs to the difference- 
equation filter description (equation (58)). A sample result of apply- 
ing a pseudo-random noise sequence to the present filter is given as 
figure 6. This figure indicates the types of "random'' time histories 
that can be developed from filters of the form (or similair form) used. 

CONCLUDING REMARKS 

The objective of this thesis has been to present an outline of the 
methods or techniques which have become available to represent crew 
motion disturbances in a manned spacecraft. The wide varieties of 
payloads, crew sizes, and control systems contemplated for future space 
missions require that updated methods of control system sizing, design, 
and analysis be utilized. The techniques included herein highlight both 
ane^.ytical and experimental approaches. Of perhaps greatest utility of 
these techniques is the z-transform application, developed here specifi- 
cally for time-domain digital computer applications requiring a variety 
of possible crew motion disturbance inputs. An example illustrates the 
ease of implementation of this method, which ease should enhance its 
utility in any future studies. 

Several other crew motion disturbance modeling techniques have been 
presented. The three deterministic methods discussed have been applied 
to ground-based simulations by a few investigators; in addition, the 
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dynamic model of man has been investigated in a flight experiment. Of 
the thre. so-called stochastic methods presented, the filtered reuidom 
noise technique has been perhaps the most useful, in frequency domain 
studies of the Skylab experiment pointing control system. Some time 
domain applications of this technique, in conjunction with pure Fourier 
transform methods have also been utilized. 

An appendix is included in this thesis; it derives a set of basic 
equations of motion for a spacecraft with applied distxirbance which, 
although requiring numerical solutions, can be used with many of the 
crew motion disturbance models (as well as with other types of defined 
disturbances ) . 

It is anticipated that results of the reseeurch efforts reported 
herein would be most useful in initial sizing studies of attitude con- 
trol systems; of almost equal utility would be their incorporation into 
existing simulation programs (both all computer and computer/hardware 
hybrid). Finally, it should be noted that the results contained in 
this thesis are not intended to be all-inclusive; no tables of para- 
meters to use with various models are given, nor are specific inplemen- 
tation schemes suggested. It is felt that a more diverse treatment of 
the subject of crew motion disturbMce modeling has resulted, with 
enough detail to facilitate applications to a variety of control system 
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Appendix 


Development of Spacecraft Dynamics Equations 


In an investigation of astronaut crew motion disturbance effects, 
relatione between disturbemee force and moment inputs and spacecraft 
response outputs are needed. Spacecraft attitude enn be eiqpressed 
in terms of Euler angles or direction cosines which in turn are 
expressible as integral functions of angular velocities about space- 
craft body- fixed axes. In order to obtedn attitude response 
information, then, it is first necessary to determine these angular 
velocities, or body rates. This appendix presents a derivation of 
the spacecraft dynamic equations, with one arbitrary moving mass, 
based on momentum considerations (refs. 23 and 2h)» A vector approach 
will be used for its economy in notation. Figure 7 illustrates the 
vector relationship between the spacecraft body-fixed axis ^ ^ 

and ^.iertial axes Xj Z^. 

The anguleu: momentum of a spacecreift, defined in coordinates 
of the body-fixed axis system (origin at 0 in figure 7 ) is the 
smn of the angular momenta of all of its mass elements, or 


H * ^ r , X m, R , 
o mj J mj 


(Al) 


and the momentum of its mass center ia 


where 




m. 


r X m R 
ms s ms 


(A2) 
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The total .ingular momentum about the spacecraft mass center is then 


H » H - H 
cm o s 


(A3) 


or 


H »\ ^ r . X n, (R + r . ♦ w x r , )- 
cm mj j o mj mj 


(Al») 


r X m (R + r + wXr ) 
ms s o ms ms 


where 5 is the angular velocity of the axis system. 

The coor-'lnates of the spacecraft mass center, r , are given by 

ms 


= ^“jLlaL = ^a.1 


m, r 


ms 


53m, 


m 


(A5) 


Equation (aU) can be rewritten usint the relation (AS) as 


H ■N"' r .n,xR - m R . m. i 

cm mJ J o “ — ’'so mJ J 


- r 


mJ ms 


m r 3 T 4 >‘n, (oi*r ,) - "r x m (wxr ) 

a ms mJ J mj mr s ms 


(a6) 


or 


H ^ ^r ,xm. r,~7 Xm r r . x»« (Cxr ,) - r x 

cm mJ J mJ ms s ms mj J mj ms 


®- (wxr ) 

8 SUB 


(AT) 


For convenience, first consiaer the spacecraft as a rigid body; r 

r ■0, and a rigid-body angular momentum can be expressed using 

m8 

equation (AT) as 


mJ 
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r ,X m, {(u X r . ) - r x m (uxr ) 
inj J nj ms s ms 


or by using elementary vector analysis 


■]C‘J [‘'■-o • V)‘“ ■ ^ V • -“>'J 


Equc*ion (A9) can be expanded into Cartesian coordinates, by 
r cognizing certain quantities defined as moments of inertia, the 
rigid-body momentum can be written in matrix form as 


I -I -I 
X xy xz 


-I I -I 
xy y yz 


-I -I I 
xz yz z 


(AlO) 


It should be noted that the moments of inertia are about axes 
through the center of mass of the spacecraft. 

Consider the origin 0 to be at the center of mass 
of the spacecraft (taken without any moving masses); then the 
consideration of some number of moving masses, k, permits the 
expression of (say) the x centroidal axis moment of inertia as 




jsn-k+l 


(All) 



Here the first term on the right-hand side of equation (All) 
represents the spacecraft moment of inertia, exclusive of any moving 
masses, and can be assumed constant for many analyses. The second 
term is the moment of inertia due to H moving masses, referred to the 
origin 0. The third term is the inertia shift to the overall 
mass center of the spacecraft and moving masses (crew, for e xampl e) . 

For convenience in subsequent development, assiime only one 
moving mass m located at a distance r from the origin 0. The 
moment of inertia about the x-axis (centroidal) can nov be written 
as 

^x “ ^x.o “s ^^s * 

where I is the first term on the right side of equation (All). 

The distance from the origin 0 to the mass center, is now 

m r 

(A13) 

3 

and 

r. n ^ 

X X ,0 D 


Defining a mass factor Q as 


m(m - m) 
s 


m 

s 


(A15) 


the x-axis moment of inertia can be written 


= I. 


x,o 


(y' 




(Al6) 
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All moments and products of Inertia can be expressed in a 
sisiilar manner: 


^x 

= 

I,. 9(y' ♦ 

z*) 

I 

y 

X 

ly. <1(»' ♦ 

z*) 

I 

z 

= 


y*) 

I 

xy 

X 

'xy-o * 


I 

xz 

= 

I . + Q(xz) 

xz*o 


lyz 

s 

I » ■•• Q(yz) 
yz*o 



(A17) 


Again considering all mass elements, from equation (a 1) is obtained 


d H 


b = d 

dt dt 


\ ' r , X m, R , 
iT-j J mj 


(Al8) 


d H 
c 

dt 




(A19) 




(A20) 


d H 
0 

dt 


m. ‘^^mj ^ ^o +\ ' r , m R 
J Lmi Bj J aj 


(A21) 


d H 
c 

dt 


m dr R ^ is 

3 _^X o **mj 


(A22) 
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Similarly, from equation (A2) 


d H. 


s SE m 


dt 


^o + r ^ ®s '^ms 

— 


From equation (A3) 


d H 


cm 


d H - d H 


dt 


dt 


dt 


d H _ 

cm => r/m.R.-r «m 
^ j mj ms : 


R 

s ms 


Newton's force equation, through the mass center is 


F = m R 
s s ms 


Equation (A2^) can now be written 


d H 


cm 


dt 


mj J mj 


r X F 
ms s 


Further, the total moment about the origin 0 is 


M mA . 

o mJ J mJ 


Equations (A27) and (A28) can be combined to yield 


M = dH+r>^F 
0 cm ms s 

dt 


(A23) 


(A2U) 


(A25) 


(A26) 


(A27) 


(A28) 


(A29) 
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For most problems involving spacecraft rotations (of principal 
concern in this development) the external force applications linked 
vith mass centei shifts can be neglected4 Thus» the applied moment 
about 0 will equal the derivative of the angular momentum about the 
overaill mass center: 


M 

o 


d H 


cm 


dt 


(A30) 


Comibining and rearranging equations (AT) and (a 8) yields 

mj j mj ms s ms rb vA31) 


H 


cm 


and (considering only one moving mass m at distance r from 0) 
equation (A30) can be expanded as 


M = d H 
o cm 

dt 


d (rXmr-r Xm r -►H.) (A32) 

— CIS s ms rb 

at 


Using equation (A13) to substitute for and equation (A32) 


becomes 


M 


Also 


o dt 


M 


d r mr m . „ 1 

-T7 — rxmr - — xm r + H , 

dt m s ^ j 


’ ^ [' 
by making use of equation (A15)> 

Expanding equation (A3^) results in 
*Q^rXr-u)X(rXr)^+ H 


Q(r ^ r) + 


“rt, ] 


(A33) 


(A3U) 


+ w X H 
rb rb 


(A35) 
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Novj by substituting for (and by making use of equation 

(AI 2 ), equation (A35) can be expanded to yield the components of 

M about the X Y • and Z - axes: 

O ODD 


.. . • ^ 
M = Q [(yz - zy) ♦ oj (xy - yx) - w (zx - xz)] + (I w + w I - 

I u) “U) I — I u) — u)I (lu> ^xz^x 

xy y y xy xz z z xz y ^ z x * * ^ 

U) (lU) - Iw-Icii) 

2^ y y yz 2 xy x 


(A36) 


® Q [{zx - xz) + U) (yz - zy) -w (xy - yx)l + (I u ♦ w I - 

^ 2 X yyyy 

% • 

I UJ— U)I — I CO ^tol )+“U'(lC0 — I (0 — I Ul) — 

yz 2 z yz xy X x xy ' z^ x x ^xy y xz 2^ 

U)(lu) -Ico-Ico) 

X ' z 2 XZ X yz y' 


z z z z 


» Q C(xy 

- yx) + 

u (zx 

X 

• 


a 

I ( 1 ) 

w I - 

I u 

X 

K 

X xz 

yz y 

1 

3 

M 

3 

I b> ~ I U 

y X X 

V 

XZ z 

If the moments 

«x* «y- 

and M 

z 


y yz X y y 


- - I u ) - . 

yz z xy X y 


X, y, z are expressed as fvmctions of time, then eqxiations (A36) 


can he nximerlcally Integrated to find the body rates u ,u> ,and u> . 

X y z 


52 



References 


1. Langley Research Center Staff: "A Report on the Research and 

Technological Problems of Manned Rotating Spacecraft." NASA 
Technical Note D-I 30 U, August 1962, pp 71-83 

2. Kurzhals, Peter R. , wd Keckler, Claude R. : "Spin Dynamics of 
Manned Space Stations." NASA Tech. Report R-133* Deconber 1963. 

3 . Kurzhals, Peter R. , and Shearin, John G. : "An On-Off Control 

System for a Rotating Manned Spacecraft." NASA Tech. Note 
D- 2726 . March I 965 . 

4. Kurzhals, Peter R. , and Grantham, Carolyn: "A Syston for Inertial 

Experiment Pointing and Attitude Control." NASA Tech. Report 
R- 2 U 7 , August 1966 . 

5 * Martz, C. William, and Grantham, Carolyn: "A Passive Balance 

for a Class of Rotating Spacecraft. " NASA Tech. Note D-6924, 
October 1972 

6. Groom, Nelson J. , Shaughnessy, John D. , and Nene, Vilas D. : "On 
the Stability and Pointing of an Attached Double-Gimbal Experi- 
ment P€u:kage.” NASA Tech. Note D-6676, April 1972 

7- Hendricks, T. C. , and ichnson, C. H. : "Stochastic Crew Motion 

Modeling" Journal of Spacecraft, Vol. 8, No. 2, Febniary 1971. 

8. Shaxighnessy, John D. , Groom, Nelson J. , and Nene, Vilas D. : 
"Effects of Certain Control-System Nonlinearities on Stability 
and Pointing of an Attached Double-Gimbal Experiment Package 
in Presence of Random Crew-Motion Disturbances." NASA Tech. 

Note D-72U7, May 1973 

9 . Conway, Bruce A. : "Development of Skylab Experiment T-013 

Crew/Vehicle Disturbances." NASA Tech. Note D 658U, January 1972. 

10. Hanavan, Ernest P., Jr.: "A Mathematical Model of the Human 

Body." AMRL-TR-64-102 , U.S Air Force, October 1964. 

11. Fuhrmeister, W. F., and Fowler, J.L. : "Experimental Study of 

Dynamic Effects of Crew Motions in a Manned Orbital Research 
Laboratory (MORL)." NASA Contractor Report 66186, October 1966. 

12. Tewell, J. R. , and Mxurish, C. H. : "Engineering Study and 

Experiment Definition for an Apollo Applications Program Experi- 
ment on Vehicle Disturbances Due to Crew Activity." NASA Contra- 
ctor Report 66277. March I967. 


53 



13. Murriah, C. H. , and Smith, G. W. : "Apollo Applications Program 

Crew Motion Experiment - Program Definition and Design Development." 
NASA Contractor Report 66599, March 1968. 

lU. Papoulis, Athanasios: "Probability, Random Variables, and 

Stochastic Processes." McGraw-Hill Book Company, New York, 1965. 

15. Blackman, R. B. , and Tuk^, J. W. : "The Measurement of Power 

Spectra.” Dover Publications, Inc., 1959- 

16. Jury, E. I.; "Theory and Application of the z-Transform Method.” 

John Wiley and Sons, Inc., New York, I96U. 

17. Smyth, Michael P.t "Linear Engineering Systems." Pergamon Press, 
Inc. New York, 1972. 

18. Kuo, Benjamin C.: "Analysis and Synthesis of Sampled-Data Control 

Systems.” Prentic-Hall, Inc., Englewood Cliffs, N.J., 1963. 

19. Gupta, Someshwar C. , and Hasdorff, Lawrence: "Fundamentals of 

Automatic Control." John Wiley and Sons, New York, 1970. 

20. Otnes, Robert K. , and Enochson, Loren: "Digital Time Series 

Analysis." John Wiley and Sons, New York, 1972. 

21. Lindorff, David P. : "Theory of Sampled-Data Control Systems." 

John Wiley and Sons, New York, 1965. 

22. "Handbook of Mathematical Functions." U.S. Department of Commerce, 
Applied Mathematics Series, 196i». 

23. Kurzhals, Peter R. : "An Approximate Solution of the Equations of 

Motion for Arbitrary Rotating Spacecraft." NASA Tech. Report 

R-269, 1967. 

24. Conway, Bruce A. : "An Analysis of the Operation of a Mass Measur- 

ing System in an Orbiting Spacecraft." NASA Tech. Note D-5039, 
February, 1969. 


54 



